Evolution of spin entanglement and an entanglement witness in multiple-quantum 

NMR experiments 
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We investigate the evolution of entanglement in multiple-quantum (MQ) NMR experiments in 
crystals with pairs of close nuclear spins-1/2. The initial thermodynamic equilibrium state of the 
system in a strong external magnetic field evolves under the non-secular part of the dipolar Hamil- 
tonian. As a result, MQ coherences of the zeroth and plus/minus second orders appear. A simple 
condition for the emergence of entanglement is obtained. We show that the measure of the spin 
pair entanglement, concurrence, coincides qualitatively with the intensity of MQ coherences of the 
plus/minus second order and hence the entanglement can be studied with MQ NMR methods. We 
introduce an Entanglement Witness using MQ NMR coherences of the plus/minus second order. 

PACS numbers: 03.67.Mn, 03.67.-a, 75.10.Pq, 82.56.-b 

I. INTRODUCTION 



Entanglement [l[ is the key concept in Quantum Information Theory. It has played a crucial role in experiments 
on quantum computing and quantum teleportation. This resulted in intensive interest to the physics of entanglement 
f"| ■ from both theorists and experimentalists [3, H, 0, IE] ■ 

Q-f Entanglement is detected with the help of a so-called Entanglement Witness (EW). By definition, EW is an 
observable which has a positive expectation value for separable states and negative for some entangled states @. In 
particular, internal energy Q and magnetic susceptibility @ were used as EW in some cases. In this paper we propose 
a new type of an Entanglement Witness, the intensity of multiple quantum coherences in spin systems. This quantity 
^-i is accessible in NMR experiments and thus opens a new approach to probing entanglement with highly advanced 
NMR techniques. 

In the present work we focus on the simplest relevant system, a pair of spins s = 1/2 coupled by the dipole-dipole 
interaction in the conditions of the multiple-quantum (MQ) NMR experiment Here the initial thermodynamic 
density matrix describing the interaction of the spins with the strong external magnetic field is subjected to the 
irradiation by the specially tailored sequence of resonance rf-pulses. The anisotropic dipolar Hamiltonian oscillates 
i rapidly when the period of the sequence is less than the inverse dipolar frequency. The spin dynamics of the system 
is described by the averaged Hamiltonian which is responsible for the emergence of the MQ coherences of the zeroth 
^2 , and plus/minus second orde rs I l0l| . It is evident that the initial state of the system is separable. However we show 
with the Wootters criterion [Tl[ that the entangled state emerges when the intensity of the MQ coherence of order 
2 (-2) exceeds the exactly calculated threshold depending on the external magnetic field and the temperature. Thus 
the intensity of the MQ coherence of the second order, which is the observable in MQ NMR experiments, serves as 
^ | EW for spin systems. 

II. MQ DYNAMICS OF A DIPOLAR COUPLED SPIN PAIR AT LOW TEMPERATURES 

We consider a two-spin system in a strong external magnetic field Hq. The thermodynamic equilibrium density 
matrix, po; of the system is 

where u>o — ^Hq (7 is the gyromagnetic ratio), T is the temperature, I a = I\ a + l2 a , and Ij a (J — 1,2; a = x, y, z) is 
the projection of the angular spin momentum operator of spin j on the axis a, and Z is the partition function. 

The MQ NMR experiment consists of four distinct periods of time: preparation, evolution, mixing, and detection 9fl. 
MQ coherences are created by the multipulse sequence consisting of eight-pulse cycles on the preparation period 9(. 
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In the rotating reference frame [12j, the average Hamiltonian, Hmq, for the two-spin system describing the MQ 
dynamics at the preparation period can be written as Q 



H 



MQ 



= b If J 2 " 



(2) 



where b — (7 2 ^/{2rf 2 })(l — 3 cos 2 #12) is the coupling constant between spins 1 and 2, j-12 is the distance between 
spins 1 and 2, and #12 is the angle between the internuclear vector 7^2 an d the external magnetic field Ho', if and J~ 
(j = 1,2) are the rasing and lowering operators of spin j. 

The two-spin Hamiltonian Hmq can be diagonalized with the transformation (in the standard basis 
{|00),|01>,|10},|11}}) 
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and the density matrix, p(r), at the end of the preparation period is 



p(r) = e- lHM ^ T p e lHM ^ T - ' 



2(1 + cosh p) 



/cosh (3 + cos(2&r) sinh (3 i sin(26r) sinh j3 
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1 

\ — i sin(26r) sinh/3 cosh j3 — cos(26r) sinh (3/ 



(4) 



where (3 = tkOo/(kT), The diagonal part of the density matrix of Eq. (j4|), p( ) ( T )i is responsible for the MQ coherence 
of the zeroth order, and the non-diagonal parts, p( 2 ) (t), P(-2) ( r ): are responsible for the MQ coherences of the 
plus/minus second orders 0, [l(|. The intensities of the MQ coherences of the zeroth, Gq(t), and plus/minus second, 
G±2(t), orders are [13] 

G (r) = Tr (p (0) (r)pg) (r)) , G ±2 (r) = Tr ( P(3) (r)p^ 2) (r)) , (5) 
where (r) is the diagonal part of 

P M {t) = e -* H MQTJ ze lH MQ T ^ 

and P( l 2)( T )' p\-2)( t ) are tne non-diagonal parts of the density matrix p ht {T). Using Eqs. dl])-© one can find that 

G (t) = tanh ~ cos 2 (26r), G ±2 (r) = i tanh | sin 2 (26r). (7) 

It is worth to emphasize that intensities of MQ coherences are observables in MQ NMR experiments. Eq. ((7|) shows 
that the intensities of the MQ coherences of the second order, Ga(r), and the minus second order, G_2(r), are equal. 
However, in real experiment, certain errors are present and the experimental results for G2(t) and G_2(t) are not 
the same. Some of the errors can be compensated and the accuracy can be improved if one detects the sum of these 
coherences fl4| . It is also worth to notice that the accuracy of the measurement of G 2 (r) + G_ 2 ( T ) is higher than for 
Go(t) [3]. Below we will use the sum of the MQ coherences of the plus/minus second order in order to introduce the 
entanglement witness. 



III. CONCURRENCE AND ENTANGLEMENT WITNESS IN MQ NMR EXPERIMENTS 

The initial state of the system determined by Eq. JT]) is separable. Entanglement appears in the course of the 
preparation period of the MQ NMR experiment when the MQ coherence of the second order has a sufficiently large 
intensity. In order to estimate the entanglement quantitatively we apply the Wootters criterion [ill ] . According to [llj ] , 
one needs to construct the spin-flip density matrix 



p(r) = (a y <j y )p*{T)(a y <g) <Jy) 



(8) 
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where the asterisk denotes complex conjugation in the standard basis j|00), |01) , 1 10) , 111)} and the Pauli matrix 



ly = 2I y . The concurrence of the two-spin system with the density matrix p(r) is equal to 11 1 



C = ?7iair{0, 2A — Ai — A2 — A3 — A4}, A = max{\\, A2, A3, A4} (9) 

where Ai, A2, A3, and A4 are the square roots of the eigenvalues of the product p(r)p(r). Using Eqs. JU), (0, © one 
obtains 



Ai. 2 ; —p , A 3 , 4 - - 2 g • (10) 



'1 + sin^(2&r) sinh^ f3 ± | sin(2&r)| sinh/3 
As a result, the concurrence, G, is 



4 cosli § 4 coslr 2 



c _ |sin(2&r)|sinh/3-l 
2 cosh 2 f 

The entangled state can appear only at sinh/3 > 1 when the intensity of the MQ coherence of the second order has 
the maximal value. This condition means that the entanglement appears at temperatures 

t< *a* . (12) 

fcln(l + \/2) 

If one takes ojq = 27r500 • 10 6 s _1 the entangled state emerges at the temperature Tg « TJmK. It is interesting to 
notice that in a linear chain of dipolar coupled nuclear spins in the thermodynamic equilibrium state, entanglement 
appears only at microkelvin temperatures [la ]. 

The simple connection between the concurrence, G, and the intensities of the MQ coherences of the plus/minus 
second orders, G±2 (•?"), can be found from Eqs. (|7|), ©: 



C = ytanh |[G 2 (r) + G_ 2 (r)] - 2 J^ & - (13) 



Thus, entanglement is possible only when 



G 2 (r) + G_ 2 (r) > (14) 
2smh/3cosli § 

and Entanglement Witness (EW) can be introduced as the following 

EW = 2^L>?r mT)+GMT}} - (15) 

In the initial moment of time G2(0) + G_2(0) = 0, EW > and the considered system is in a separable state. In the 
course of the MQ NMR experiment the intensity of the MQ coherence of the second order grows and the entanglement 
witness, EW, changes its sign. It means that an entangled state appears. According to Eq. ([7]) the intensities of the 
MQ coherences periodically change in time. The sign of EW changes also periodically. Thus separable and entangled 
states change periodically in the considered system. The time evolutions of the MQ coherences of the zeroth and 
second orders together with the corresponding concurrence are represented in Fig. Q] at (3 = 3. One can see that 
the concurrence is close to the sum of the MQ coherences of the plus/minus second orders, G2(t) + G_2(r), at 
almost all durations of the preparation period of the MQ NMR experiment. At large j3 (small temperatures) the 
expression [2 sinh ft cosh 2 tends to zero and the maximal value of G2(r) + G_2(t) tends to one. This means 
that the concurrence coincides with the maximal value of G 2 (r) + G_2(t) at small temperatures. The corresponding 
dependencies of the concurrence and the maximal value of the sum G2(r) + G_2(t) on (3 are given in Fig. [2l We can 
conclude that the entangled states appear in MQ NMR experiments at sufficiently small temperatures. In contrast 
to Ref. Q we study entanglement in a system of nuclear spins(not electron ones). Such systems are more robust 
to decoherence which causes the loss of the quantum information which was obtained during quantum computation. 
A problem, related to ours, was studied in Ref. [1 €31 ] - That work focuses on the high temperature reg ime in which 
entanglement is absent. The prediction of the existence of an entangled state at high temperatures [161 ] is an artifact 
of an incorrect choice of the initial density matrix. 
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FIG. 1: The dependencies of the MQ coherences and the concurrence on the time of the preparation period, r, of the MQ NMR 
experiment at /3 = 3. The coupling constant is equal to b = 27rl307s~ 1 ; solid line - concurrence; dashed line - intensity of the 
MQ coherence of the zeroth order; dash-point line - (?2(t) + G_2(t) (see the text). 
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FIG. 2: The dependence of the concurrence (dashed line) and the maximal value of G'2(t) + G_2(t) (dash-point line) on the 
parameter j3. The solid line describes the function [2 sinh (3 cosh 2 f Here the coupling constant is equal to b = 27rl307 s -1 The 
entangled state emerges at temperatures less than Te = huo/(kf3E). 



IV. CONCLUSION 



MQ NMR experiments can be used for the analysis of entangled states in spin systems. We have introduced an 
Entanglement Witness using observable intensities of the MQ NMR coherences of the plus/minus second order and 
analyzed entanglement in term of the Wootters criterion. Entangled states emerge when the sum of intensities of the 
MQ coherences of the plus/minus second order exceeds an exactly calculated threshold depending on the external 
magnetic field and the temperature. MQ NMR can be considered as a new method for obtaining entangled states in 
spin systems. 
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